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ABSTRACT 



Context. Mass loss from massive stars forms an important aspect of the evolution of massive stars, as well as for the enrichment of 
the surrounding interstellar medium. 

Aims. Our goal is to predict accurate mass-loss rates and terminal wind velocities. These quantities can be compared to empirical 
values, thereby testing radiation-driven wind models. One specific topical issue is that of the so-called "weak-wind problem", where 
empirically derived mass-loss rates and (modified) wind momenta fall orders of magnitude short of predicted values. 
Methods. We employ an established Monte Carlo model and a recently suggested new line acceleration foiTnalism to solve the wind 
dynamics more consistently. 

Results. We provide a new grid of mass-loss rates and terminal wind velocities of O-type stars, and compare the values to empirical 
results. Our models fail to provide mass-loss rates for main-sequence stars below a luminosity of log(L/Lo) = 5.2, where we appear 
to run into a fundamental limit. At luminosities below this critical value there is insufficient momentum transferred to the wind in the 
region below the sonic point in order to kick-start the acceleration of the flow. This problem occurs at almost the exact location of the 
onset of the weak-wind problem. For O dwarfs, the boundary between being able to start a wind, and failing to do so, is at spectral 
type 06/06.5. The direct cause of this failure for 06.5 stars is a combination of the lower luminosity and a lack of Fe v lines at the 
base of the wind. This might indicate that - in addition to radiation pressure - another mechanism is required to provide the necessary 
driving to initiate the wind acceleration. 

Conclusions. For stars more luminous than 10^^ Lq, our new mass-loss rates are in excellent agreement with the mass-loss prescription 
by Vink et al. 2000 using our terminal wind velocities as input to this recipe. This implies that the main assumption entering the method 
of the Vink et al. prescriptions - i.e. that the momentum equation is not explicitly solved for - does not compromise the reliability of 
the Vink et al. results for this part of parameter space. Finally, our new models predict terminal velocities that are typically 35 and 45 
percent larger than observed values. Such over-predictions are similar to those from (modified) CAK-theory. 

Key words. Stars: early-type, massive, mass-loss, winds, outflows. Radiative transfer 



1. Introduction 

I In this article, we present predictions for mass-loss rates and ve- 
. locity structures for a grid of O-type stars, using two distinct 
' methods for solving the wind dynamics. 

Mass loss forms an integral aspect characterizing massive O- 
' type stars. Because of their short lifetimes, massive stars are im- 
' portant tracers of star formation in galaxies. Furthermore, they 
enrich the interstellar medium with metals, both during their 
lives via stellar winds, as well as when they explode at the 
very end of their evolution. In order to build an evolutionary 
framework for massive stars, it is essential to map the mass-loss 
processes during the various evolutionary stages, as the exact 
rates of mass loss greatly influence the evolutionary tracks (e.g. 
iMaededfigsTtlChiosi & MaedeiilT986l) . The effects of mass loss 
on the evolutionary tracks are at least two-fold: first and fore- 
most the stellar mass is reduced, and secondly, the rotational ve- 
locity is strongly affected, as the mass also carries away angular 
momentum (e.g. Langer 1998; Maeder & Meynet 2000). 

For the continuous stellar winds from massive stars, the 
outflow is thought to be driven by the transfer of energy 



and momentum from the radiation field to the atmosphere 
through the absorption of photons in atomic transitions. The 
exact amount of momentum and energy transfer has been 
the subject of both theoretical and obser vational stud i es for 
many decades (Lucy & Solomon 1970", 



Pauldrach et all [T986t .Puis et al., ,1996: 



Castor et al 



Vink et alJll999t iKrticka & Kubatll2004t 



ide Koter et al 



Mokiem et al 



1975 



19971 : 
I2OO7I) . 



For luminous O-ty pe stars, w ith logiL/Lo) > 5.2, the theoret- 
ical predictions of I Vink etal. ( 2000) seem to be in reasonable 
agreement with empirical mass-loss rates provided that O-stars 
are only subject to modest amounts of wind clumping (with 
clump filling factors of only 5-10). However, for objects with 
luminosities log(L/Lo) below approximately 5.2, a severe drop 
- by a factor of ~ 100 - in the empirically determined modi- 
fied wind momentum (basically a multiplication of the mass- 
loss rate and the terminal velocity) has been revealed. This prob- 
lem ha s in literature been referred to as "the weak-wind prob- 
lem" dPuls et all 1 1 9961: iMartins et al] l2005bl: IPuIs et al J 120081: 
iMarcolino et al.ll2009l) . 
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It deserves proper investigation simply because of the enor- 
mity of the effect. It is particularly important to find out whether 
the problem is caused by the mass-loss diagnostics or the pre- 
dictions, as both are also applied to more luminous stars, where 
agreement between diagnostics and predictions has seemingly 
been achieved. But how certain can we be that this agreement is 
not a coincidence if we are aware of severe problems at lower 
luminosity? 

Furthermore, we note that the oft-used mass-loss predictions 
of lVink et al.l (l2000l) are semi-empirical, in the sense that empir- 
ical values for the wind velocity structure and terminal velocity 
are used as input to the modelling. In order to trust our over- 
all knowledge of the mass-loss rates from O-type stars - at all 
masses and lumi nosities - it is pivotal to further scrutinize the 
IVink et all ( l2000h assumptions, most notably that of the adopted 
wind dynamics.^ 

Recently, iMuller & Vin3 (l2008h suggested a new 
parametrization of the line acceleration, expressing it as a 
fu nction of rad i us ra ther than of the velocity gradient, as 
in ICastor et al] ( 119751: henceforth CAK) theory. The imple- 
mentation of this new formalism allows for local dynamical 
consistency, as one can determine the energy and momentum 
transfer at each location in the wind through the use of Monte 
Carlo simulations. Although the formalism was applied with 
three independent starting conditions that showed convergence 
to the same wind parameters, it has thus far only been applied to 
one object, that of an 05 dwarf. 

For the adopted line force parameterization iMiiller & VinkI 
identify an exact solution in case the medium is isothermal. We 
expand on this result by also accounting for a temperature strati- 
fication. To allow for such a study we employ the new line accel- 
eration parameterization but solve for the wind dynamics consis- 
tently by applying a numerical method to solve for the momen- 
tum equation. 

The purpose of our study is threefold; i) to solve the wind 
dynamics nurnerically , and compare the results to those of 
iMiiller & VinkI (l2008h . ii) to compute a larger grid of dynam- 
ically derived O-star mass-loss rates and wind terminal veloc- 
ities, and d eterm ine the accuracy of the predictions made by 
IVink et all (l2000l) . and Hi) to utilize the grid in order to inves- 
tigate the weak-wind problem. 

Our paper is organized as follows. In Sect. |2] we start off 
describing the core of our method and the different methods to 
treat the wind equation. The results are presented in Sect. [3] and 
discussed in the Sect.|4| We end with the conclusions (Sect.|5]l. 

2. Method 



The method of Ide Koter et all (Il997h and IVink et al.l (1 19991) . 
applied to deri ve the mass-loss rates of O and early- 
B type stars dVink et al.l l2000i 



200l . Luminous Blue 



Variable stars (IVink & de Kotej l2002h and Wolf-Rayet stars 



dVink&d e Koter '20051). is an extension of a treatment devel- 



oped bv [Abbott & Lucy! (Il985l) . It is based on an iteration cycle 



de Koter et al 



de Koter et al 



between the stellar atmosphere model isa-wind 
1993 ) and a Monte Carlo simulation, mc-wind 
19971) . in which the energy per unit time AL that is extracted 
from the radiation field in interactions of photons with the gas, 
is computed. From this a mass-loss rate M is computed on the 
basis of which a new isa-wind model is constructed. The pre- 
dicted mass loss is the one for which the input mass-loss rate of 
ISA-WIND equals the mass-loss rate computed by mc-wind. 

As is consistently pointed out in the papers referred to above, 
the method does not address the equation of motion but uses a 



prescribed trans-sonic velocity structure. This implies that al- 
though in a global sense the method fulfills the constraint of en- 
ergy conservation, it need not hold that the actual local forces 
acting on the gas are consistent with t he for ce implied by the 
adopted velocity law. IMiiller & VinkI (|2008|) relax on this as- 
sumption and present an improved treatment of the problem in- 
troducing a new way to parametrize the line force. We first dis- 
cuss an approach presented by these authors, which we refer to 
as the "best-yS method, as it allows to link to empirically derived 
estimates of the steepness of the velocity law (characterized by 
a parameter see below). In a second step, we present solutions 
that numerically solve the wind dynamics. 

We first briefly introduce isa-wind in Sect. 12.11 emphasiz- 
ing the treatment of the heuristic velocity law, and mc-wind 
in Sect. 12.2! focusing on the determination of the mass-loss 
rate using the global energy argument. In Sect. 12.3! we reca- 
pit ulate the essent ials of the parametrization of the line force 
bv IMuller &"V^ and in Sect. 12.4! the principle of their best-/? 
method. In the following subsection we introduce our hydrody- 
namical method. Finally, Sect. I2.6l is devoted to a discussion of 
the physical conditions at the sonic point. 

2. 1. The model atmosphere 

The code isa-wind computes the structure, radiation field and 
ionization/excitation state of the gas of an outflowing stellar at- 
mosphere in non local thermodynamic equilibrium (non-LTE), 
assuming radiative equilibrium. No artificial separation between 
the photosphere and wind is assumed. The temperature struc- 
ture is treated somewhat simplified in that it results from ini- 
tial LTE based Rosseland opacities (i.e. grey). The fact that the 
temperature structure is not affected by possible departures from 
the populations from their LTE state implies that the effect of 
line blanketing is not treated self-consistently, although non-LTE 
line blocking is taken into account. Radiation transfer in spec tral 
lines is treated in the Sobolev approximation (ISobolevlll960l) . 

The input stellar parameters are the luminosity L, the effec- 
tive temperature T^s (specifying the radius R), the mass M and 
chemical abundances. The wind is described by the mass-loss 
rate M and a velocity structure, which are connected through the 
equation of mass continuity 



M - Anr v{r)p(r), 



(1) 



where p(r) is the mass density and v{r) is the velocity at radius 
r. Outside the photosphere, the velocity structure is assumed to 
follow a j6-law, i.e. 



v(r) -Voo\l- — 



(2) 



The free parameter /? is a measure of the velocity gradient. A 
low value of (3 implies that the velocity approaches the terminal 
flow velocity relatively close to the star; for a large value this 
happens only further out in the wind. The y6-law does not hold 
in the photosphere since the line force is not the dominant term 
in the equation of motion, but gravity and the acceleration due 
to the gas pressure gradient also contribute to the flow structure. 
The radius r' is a smoothing parameter that is used to connect 
the j6-law to the (quasi) hydrostatic photosphere and must as- 
sure that v{r) and its spatial derivative are continuous at the point 
where one couples the photospheric velocity law to the y6-law. 
The velocity structure in the photosphere is determined by solv- 
ing the non-isothermal equation of motion, neglecting line radi- 
ation pressure and assuming that continuum radiation pressure 
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is the result from Thomson scattering only. An inner boundary 
velocity (or density) is chosen, which may be used to tune the 
total Rosseland optical depth of the photosphere and wind (see 
below). 

The wind is assumed to be homogeneous, i.e . the outflowing 
gas is not clumped (but see iMuiires et al.ll201lT) . and the termi- 
nal velocity is chosen to be 2.6 times the effective escape ve- 
locity from the stellar photosphere, which is in reasonable con- 
corda nce with empirically determined terminal ve locities of O- 
stars dLamers et alJll995UKudritzki & Pulsll2000() . The base of 
the photosphere is positioned at a Rosseland optical depth of 
about 20-25 and the wind extends out to 20 



2.2. The Monte Carlo method mc-wind 

The code mc-wind uses the model atmosphere structure com- 
puted by ISA-wind to determine the total amount of energy that is 
transferred from the radiation field to the wind - in interactions 
of photons with ions in the gas - by means of a Monte Carlo sim- 
ulation of the trajectories of photon packets emitted at the base of 
the photosphere and escaping through the outer boundary of the 
model. Each photon can travel an optical depth weighted (ran- 
dom) distance to a point of interaction. This point is determined 
by taking into account all the opacity the photon encounters on 
its path, so it includes contributions from both lines and con- 
tinua. At the point of interaction the type of interaction is deter- 
mined, using proper weighing functions (IVinketal.1119991) . The 
possible interactions are thermal absorption and emission, elec- 
tron scattering and line scattering. The interaction is assumed 
to be coherent in the co-moving frame of the ion. In the ob- 
servers frame, however, energy can be exchanged from the ra- 
diation field to the gas (or vice versa). It is traced which ion is 
involved in the interaction, such that, for instance, the contribu- 
tions to the radiative force can be dissected and identified. This 
provides us with a powerful tool to study the nature of the line 

force at each location in the wind. 

T he radiative force per unit mass equals jAbbott & Lucvl 
Il985h : 



90000 



l^rad 



I dL 
M dr' 



(3) 



where dL is the amount of energy lost by the radiation field in a 
layer of thickness dr. 

Once the total amount of energy transferred to the wind is 
known, the mass-loss rate that can be driven/or the density and 
velocity structure of the adopted isa-wind model can be calcu- 
lated. Neglecting enthalpy: 



(4) 



where AL is the total amount of energy lost by the radiation field 
and 



2GM. 



(5) 



is the Newtonian escape velocity from the stellar surface. G is the 
gravitational constant. A new isa-wind atmosphere, adopting the 
mass-loss rate as determined in mc-wind, is computed followed 
by a new Monte Carlo simulation. This procedure is repeated 
until the input mass-loss rate of mc-wind equals the output mass- 
loss rate. Although the mass-loss rate that is predicted in this 
way reflects that in a global sense the energy that is needed to 
drive the wind is indeed extracted from the radiation field, it does 
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Fig. 1. The line force (blue crosses) as predicted by mc-wind in 
the first iteration step. A fit (black dotted line) using Eq. |6] to 
represent the line force is overplotted. Note the modest scatter 
on the Monte Carlo results due to noise. 



not mean that the input line force (implied by the velocity law) 
equals the output line force from the Monte Carlo simulation 
locally, i.e. the equation of motion of the wind is not solved. 

Here we improve on this situation using two methods. Both 
methods A and B require a parametrization of the line force pre- 
dicted by MC-wiND. We therefore first discuss this aspect. 

2.3. Line force parametrization 

Figure [T] shows the Monte-Carlo line force (blue crosses) as 
is produced in the first iteration step of a typical 03V star 
(L = lO^-*^ Lo, reff= 44,600K and M = 58 Mq). The Monte 
Carlo line force is determined in a statistical way and shows 
scatter Given the delicate nature of the equation of motion it 
can not be used as such and must be represented by an appro- 
priate analytical fit function. We adopt a parametrization of the 
line force as a function of radius, r ather than of optical depth, 
as opted for by ICastor et all d 1 9751) . In Sect. |42] we show that 
this leads to a more accurate numerical representation of the line 
force, at least for th e type of stars studied here. In doing so, we 
follow lMiiller & Vink (2008). who motivate 



yline 
Tad 





go{l-r,lr)yir'' 



if 
if 



r < ra 
r>ro 



(6) 



where go, r^ , and y are fit parameters to the Monte Carlo line 
force. This choice of the fit function, i.e. without any explicit de- 
pendency of the line force on the velocity gradient, implies that 
in our models the critical point is the sonic point. Figure[T]shows 
a typical result for this fit (black dotted curve). The deviations 
are (as mentioned) due to scatter in the simulation. 



2.4. Method A: Best-/3 solution 

In this section, we use the line force representation Eq.|6]to de- 
termine - after making certain assumptions - an analytical solu- 
tion of the velocity la w in the outer part of the wind, following 
iMiiUer & VinkI (l2008h . This solution can be compared to the /?- 
law (Eq. IS and used to derive u^o and the /3 value that is most 
representative. This is useful in comparing to the often applied 
;6-law. 
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We aim to find a solution of the equation of motion 



In the limit r — > oo we find that: 



du 
'} — 
dr 



2r2 



j^line 
5rad 



]_dp_ 

p dr ' 



where p is the gas pressure and 



fescN Vl - r, 



(7) 



(8) 



is the effective surface escape velocity of the star F is the con- 
tinuum radiation pressure in units of the Newtonian gravitational 
acceleration. Sufficiently far from the photosphere this term is 
dominated by radiation pressure on free electrons, i.e. r = re, 
where Fg is essentially constant for early-type stars. Close to or 
in the photosphere, the acceleration due to free-free, bound-free 
and bound-bound processes may compete with electron scatter- 
ing and should, in principle, be considered in Eq.[8] For our best- 
jS solution, however, we assume a constant continuum accelera- 
tion, which we set to Fg. Substituting the equation of state for an 
ideal gas and using Eq. [T] the term {l/p)dp/dr can be written 
as: 



1 dp 
p dr 



dv 
V dr 



2a~ 
r 



k dT 
m dr ' 



(9) 



where k is Boltzmann's constant, m the mean particle mass and 
a{r) is the local sound speed, given by: 



a - 



(10) 



We assume the wind to be isothermal, such that the sound speed 
is constant. The equation of motion can now be rewritten as 



Iflo 



dv 
dr 



r 



' rad ' 



(11) 



where Uo is the isothermal sound speed at the effective temper- 
ature of the star Equation [TT] is a critical point equation, where 
the left- and right-hand side vanish at the point v{rs) = Uo, i.e. 
where is t he radius of the sonic point. It yields several types 
of solutions. iMiiller & Vin3 (l2008h show that for the isothermal 
case and a line force as described in Eq.|6l analytical expressions 
for all types of solutions of Eq.[TT|can be constructed by means 
of the Lambert W function (as for a further discussion of the so- 
lution of^q. [TT] containing an additional centrifugal term, see 
iMiillerilMll) . Even for the interesting trans-sonic case of a stel- 
lar wind, the analytical solution has an intricate shape. However, 
a useful approximate wind solution for the velocity law can be 
constructed if the pressure related terms 2a^/r and a/v can be 
neglected. We note, however, that at the sonic poi nt the contribu- 
tion o f the two pressure terms is non-negligible dMiiller & VinkI 
|2008|) . After some manipulation one finds that the approximate 
velocity law is given by: 



^ r ra(l+y)\ rl 



(12) 



where C is an integration constant. From this equation, the ter- 
minal wind velocity can be derived if the integration constant C 
can be determined. This can be done assuming that at radius r^ 
the velocity approaches zero. This yields 



C = - 



r^ 



(13) 



2 

esc 



(1 +r) 



(14) 



The terminal velocity Voo can also be determined from the 
equation of motion. At the critical point, the left-hand and right- 
hand side of Eq. [TT]both equal zero. Introducing in relation 
to go as expressed in Eq. [14] we find 



/ r, y r, (Vesc , \ 



y2 

"^esc 



(15) 



A direct comparison to the y6-law can be made for the super- 
sonic regime of the wind and results in 



1 + y 



(16) 



Given the assumptions made in this derivation, this result is only 
approximately correct. 

The procedure that is followed to obtain the best-/? solution 
is that in each Monte Carlo simulation the values of go, r^, and y 
are determined by fitting the output line force. Using these values 
and the current value of the sonic pointradius, Eqs. [T4l[T5] and[T6l 
are used to determine Uoo and Doo derived from Eq.[T5] the mass 
loss predicted in mc-wind, and the expression derived for/3 serve 
as input for a new isa-wind model. The two codes are iterated 
until converge nce is achieved. 

Following IMiiller & Vinl (|2008|) . we assume that conver- 
gence is achieved when the values for Wco derived from Eqs. [14] 
and [15] agree within 10 percent. This implies that our predicted 
terminal velocities have at least this uncertainty. Once the ve- 
locity convergence criterion is fulfilled all fit parameters and the 
values for M and the sonic point radius will be stable to within 
five to ten percent. 

2.5. Method B: Hydrodynamic solution 

The accuracy of the best-/? solution hinges on the assumptions 
that the wind is isothermal and that the Eddington-factor F is 
constant (here taken to be equal to Fg). It may be expected that 
these assumptions have an impact on the velocity structure near 
the sonic point, which is where the mass-loss rate is set. To as- 
sess this impact and to improve on the physical treatment, we 
devise a numerical solution of the equation of motion (Eq. [TTj) 
throughout the entire photosphere and wind, referred to as the 
hydrodynamic solution. 

To this end we start our solution at the critical point v - a 
and proceed both down-stream and up-stream using a 4'*^ order 
Runge Kutta method with adaptive stepsize control (Press et alJ 
1992 ). Applying I'HopitaFs rule (see e.g. Lamers & Cassinejlj 
1999h an expression can be devised to determine dv/dr at u{r^) = 
a. In order to determine the location of the sonic point we 
require 



2r? 



line 



= 0. 



(17) 



The above equation is solved numerically. 

So far, the hydrodynamic solution assumes an iso-thermal 
medium. At the sonic point the temperature gradient is very 
small, therefore the location of r^ can be reliably determined 
using Eq. [17] even if dT/dr would be taken into account. The 
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Fig. 2. A comparison between the best-y6 and hydrodynamic velocity laws for an 03 V star. Here, the best-y6 velocity law is the one 
determined by using the fit-parameters of the hydrodynamic solution. The label 'approx' in the left panel refers to the approximate 
velocity law, as given by Eq. [12] Three different radial regimes are plotted: full radial range (left panel); the region around the 
sonic point (middle panel) in which the location of the sonic point in the hydrodynamic solution is indicated with an arrow, and 
the photospheric region (right panel). Note that at high velocities both methods, as well as the approximate velocity law, yield very 
similar velocity stratifications. Near the sonic point the best-/? velocity law is steeper than the hydrodynamic velocity law. Its sonic 
point is closer to the star Given this difference, the right panel shows that in the photosphere where hydrostatic equilibrium controls 
the equation of motion, the shape of both velocity profiles is very similar. 



neglect of the temperature gradient in the hydrodynamic solu- 
tion in the region below the sonic point has a significant impact 
on the structure - for instance on the total (Rosseland) optical 
depth from the inner boundary to the sonic point. To solve this 
problem, we account for the temperature structure inward of the 
critical point. This requires an iterative procedure between the 
solution of the non-isothermal equation of motion 

jv a\dv 2a- k dT 

and the temperature structure (see Sect. 12.1b . After starting nu- 
merical integration of the velocity structure at the sonic point 
now determined by applying Eq.[T71 but using the local value of 
the temperature at the sonic point) in the down-stream direction, 
we include the dT/dr term in Eq.[T8] This implies that the loca- 
tion of the sonic point is not affected. In the up-stream direction 
the temperature gradient is negligible, and is ignored. 

Figure |2] compares the best-/? and hydrodynamic velocity 
laws for an 03 main sequence star It shows that the best-j6 so- 
lution behaves very similar to the numerical velocity law. Note 
that if one zooms in on the location of the sonic point, one sees 
that for the best-y6 velocity structure is positioned slightly 
more inward, or, alternatively, that the velocity law is steeper 
in the lower part of the wind. In the best-yS method, the abso- 
lute scaling of the velocity structure in the photosphere is based 
on the adopted velocity at the inner boundary of the model (see 
Sect. 12. lb . Therefore, only the position of as predicted by the 
hydrodynamical method is physically meaningful, albeit in the 
context of the assumption F = Fg. 

Once this iterative procedure has converged, and the non- 
LTE state of the gas is computed throughout the atmosphere, 
we iterate between isa-wind and mc-wind in the same manner as 
described in Sect. 12.41 Again, save for the fit parameters con- 
verge on an accuracy of better than ten percent in a few iteration 
cycles. For Voa we are forced to adopt an accuracy of 20 percent. 
Our predicted terminal velocities have at least this uncertainty. 

2.5.1. Remaining assumptions and uncertainties 

In the hydrodynamic solution the contribution of bound-free and 
free-free opacity to the continuum radiation pressure is ignored 
(see Sect. 12. lb . In the photosphere, the contribution of these pro- 



cesses to F is not negligible and may in fact be of the order of 

Fe. 

We use the Sobolev approximation for line radiation transfer. 
The Sobolev approximation becomes ill-founded for small ve- 
locity gradients dv/dr or velocities lower than the sound speed. 
iPauldrach etal] ( |I986|) showed that in the photosphere (where 
the velocity is very small) the line force is underestimated in the 
Sobolev approximation. Further out, in the region of the sonic 
point, the line optical depth is overestimated compared to co- 
moving frame values, implying an overestimate of the line force 
in this region and therefore an overestimate of the mass-loss rate. 

In addition to the above two sources of uncertainty to the 
balance of forces at and below the sonic point is the quality of 
the fitting function Eq. |6] in this part of the wind, that may be 
uncertain by up to a factor of two. This is not expected to be a 
big problem deep in the photosphere, as both the fit function as 
well as the simulated Monte Carlo line force are small compared 
to the radiative force on free electrons, but at the sonic point it 
might play a role. 

2.6. The line force at the sonic point: a test for the validity of 
the best-p method 

The critical point of the equation of motion is the sonic point. 
A dissimilarity between the sonic point and the critical point 
may occur when the line force is represented by an explicit 
func tion of du/dr, such as in CAK and modified-CAK the- 
ory (IPauldrach et al.|[r986l) . Although these descriptions provide 
extremely valuable insights, they do make assumptions regard- 
ing the behavior of (See Sect. 14.2b . The same applies for 
our method. Here we want to point out that Eq. [17] implies that 
- whatever the description of the line force - at the sonic point 
•?rad ~ Seff^ the prcssurc gradient terms 2a^/r and (k/m) dT/dr 
are small compared to the line force. This is characteristic for 
monotonic winds of early-type stars (see|Feldmeier & Shlosma^ 
120001 for non-monotonic flows). Here g^g - GMt(l - F)/r^ = 
Rtvl^^/2r~. This implies that for the velocity structure to be a 
physical solution it m ust be that at the son ic point gj.™ /getf - 1, 
as pointed out by e.g. ICastor et al.l (Il975h . We require from our 
best-yS solutions, that this criterion is fulfilled. If g|.™/^eff is not 
approximately equal to 1 at the sonic point, we interpret this as 
a failure of the wind to become trans-sonic due to a lack of line 
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force at the location in the wind where it is essential. Dynamical 
effects might occur, such as fall back, that are beyond the topic 
of this paper. In any case, we interpret such solutions as cases 
in which the wind cannot be initiated by line driving alone. For 
the hydrodynamical solution a failure to fulfill the above require- 
ment implies that we do not find a solution at all. 

3. Results 

3.1. Grid 

In order to study our predictions of the wind properties of O- 
type stars in a systematic manner, we define a grid of main 
sequen ce, giant and supergi ant stars using the spectral calibra- 
tion of iMartins et al.l (12005 ab adopting theoretical effective tem- 
perature scales. This calibration is based on non-LTE models 
that take into account line blanketing effects and an outflow- 
in g stellar wind. We have ap plied solar abundances as derived 
by Anders & Grevess^ ( Il989l) . consistent with the predictions of 
IVink et all (l2000l) . For all stars in the grid, we have derived the 
mass-loss rate, terminal velocity and /^-parameter The hydrody- 
namic solution does not feature a /3, rather y is the parameter that 
describes the slope of the velocity law. To better facilitate a com- 
parison between the different methods we have applied Eq.[T6]to 
convert y into a y6-value, referred to as /3y. The calculated grid 
is given in Table [T] The final column lists the mass-los s rate as 
predicted using the fitting formula of IVink et al.l (|2000|) that as- 
sume a fixed value of/? = 1, but for an input variable value of 

Vco - 2.6 Vesc- 

Figures [3] H and |5] show wind properties as a function of 
effective temperature, for dwarfs, giants and supergiants respec- 
tively. 

We present all the results of the best beta method, i.e. prior 
to applying the requirement defined in Sect. l2.6l that at the sonic 
point the acceleration due to the line force should approach the 
effective gravity. Having pointed this out, we first draw attention 
to the striking behavior in our best-yS predictions of dwarfs. In the 
direction of decreasing temperature, the terminal flow velocity 
drastically increases for spectral types 06.5 or later. For giants 
the 07 star shows a similar behavior. We argue below that this 
behavior reflects the failure of the wind to become supersonic, 
therefore we interpret these solutions to be non-physical. 



3.2. Early 0-stars (spectral types 03 through 06) 

Let us, however, first focus on stars of spectral type earlier than 
06.5. The two methods give quite comparable results. The best- 
y6 method predicts M values that are higher by up to ~0.1 to 
0.3 dex in all cases, i.e. dwarfs, giants and supergiants. The best- 
jS terminal flow velocities are ~ 10 to 20 percent lower compared 
to the hydrodynamic solutions. These differences can be under- 
stood by focusing on the velocity structures near the sonic point. 
In the best-/? solution the velocity law is steeper in the region 
near the sonic point, therefore the sonic point is closer to the 
photosphere. This leads to a higher mass-loss rate and lower ter- 
minal velocity. The absolute value of the terminal velocity and 
the ratio of u^o to the effective escape velocity as a function of 
temperature will be compared to observations in Sect. |4] Typical 
error bars on the Uoo determination are 10 percent for the best-/? 
solutions (see Sect. \2.4\ and 20 percent for the exact solutions 
(see Sect. 12.5b due to Monte Carlo noise on the line force (see 
also Fig.[TJ. 

The slope of the velocity law in the best-/? solution increases 
slightly with luminosity class, from typically 0.85 in dwarfs, to 



0.95 in giants, to 1.0 in supergiants. In the hydrodynamic models 
(method B) the /3y value is typically 0.05-0.10 lower than the 
corresponding best-yS solution. 



3.3. Late 0-stars (spectral types 06.5 through 09.5) 

Figures|3]and|4]show that for spectral type 06.5 the terminal ve- 
locity of dwarfs and giants suddenly peaks, relative to spectral 
type 06. We investigate this behavior in more detail in Fig.|6]in 
which the line force from the Monte Carlo simulation is plotted 
in the region around the sonic point for the best-/? solutions of the 
dwarf 06 and 06.5 star. For the 06 star, the line force at the base 
of the wind (below the sonic point) rises steeply. At first the dom- 
inant contributors are iron lines, notably from Fe v. The ensem- 
ble of transitions mainly occur between excited states, although 
some are from meta-stable states that are relatively strongly pop- 
ulated. Further out, the iron contribution levels out (at ~ 30%) 
and other elements start to contribute to the force, such as car- 
bon, nitrogen, sulfur, argon and nickel. The contribution of reso- 
nance lines of carbon and nitrogen at the sonic point amounts to 
~ 20%. Note that at the sonic point the §|™/geff ~ 1 condition is 
nicely fulfilled for the 06 V. 

For the 06.5 star this is not the case. Here the line force at 
the base of the wind (below the sonic point) rises only gradually. 
The difference with the 06 V star is that in this region iron is 
mainly in the form of Fe iv, which for this particular spectral flux 
distribution is less efficient in absorbing stellar flux than are Fe v 
line^. Therefore the velocity structure will be shallower, limit- 
ing the potential of other elements in contributing to the force. 
As a result the sonic point starts to shift out to larger radii, and 
we find that at the sonic point the cumulative line acceleration is 
some 40% less than the effective gravity. We therefore interpret 
this outcome as a failure of the wind to become supersonic at r^ 
and do not consider it to be a physical solution. 

The best-/? solutions where we clearly encounter this prob- 
lem have brackets placed around the predicted wind properties 
as listed in Table [T] These include all the dwarf stars of spectral 
type 06.5 or later They are to be considered non-physical. 

The supergiants do not suffer from this problem. In all cases 
the g|.™/^etf ~ 1 was reached at the sonic point and we consider 
them physical solutions. The terminal velocities for the 06.5 1 
to 09.5 1 scatter by about 20%, with a small hint that here also 
the 06.5 star has a higher v^o. The latter occurs because elements 
such as silicon, iron and sulfur add to the line force in the outer 
wind along with the normal contribution of carbon, nitrogen and 
oxygen. 

The value of yS for the late spectral types increases to 1 .05 
from 1.0 for earlier spectral types. The jSy values associated to 
the hydrodynamic solutions increase marginally compared to 
that in early-O stars. 



4. Discussion 

In discussing our results we first compare with previous theo- 
retical predictions for mass-loss rates and terminal velocities in 
sections 14711 and W2\ We compare to observations in section |43] 



' We note that a similar situation occurs at spectral type Bl, where 
the relative ly inefficient Fe rv lines are replaced by the more effective 
Fe m lines jVink et all 19991) . 
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Table 1. Model parameters, following iMartins et al.l (12005 ah . and pre dicted wind p roperties for dwarfs, giants and supergiants. 
The label "spec" indicates that spectroscopic masses are adopted from lMartins et al] Predictions give the mass-loss rate, terminal 
velocities and [} parameters for both method A (best-j6 solution) and B (hydrodynamic solution). The 11th column states whether or 
not the best-y6 solution fulfills the requirement that at the sonic point §[.™/geft - 1 (see Sect. 12.61 ). For the hydrodynamic solutions a 
failure of this requirement implies that we do not find a solution at all. The last column provides the M bv lVink et al.l (1200 lb when 
we use = 2.6 Wesc in their mass-loss recipe. The hydrodynamic solution does not provide a j3 value, but rather the fit parameter y. 
As to facilitate a comparison, we applied Eq.[T6]to convert this y into a fi. 



Model Parameters 












Method A 






Method B 






Vink et al. 


ST 


Teff loj 




logL 


R 




«csc 


logM 


Voo 




r=: 1 


logM 


Vac 




logM 




K cm s - 


Lo 


Rq 


Me 


km/sec 


log Mo/yr 


km/sec 




at (J , 


Mo/yr 


km/sec 




Mo/yr 


Dwarfs 




























3 


44616 


3.92 


5.83 


13.84 


58.34 


1054 


-5.641 


3794 


0.92 


yes 


-5.972 


4530 


0.87 


-5.375 


4 


43419 


3.92 


5.68 


12.31 


46.16 


1016 


-5.836 


3599 


0.90 


yes 


-5.929 


3973 


0.83 


-5.571 


5 


41540 


3.92 


5.51 


11.08 


37.28 


992 


-5.969 


2838 


0.84 


yes 


-6.118 


3394 


0.79 


-5.829 


5.5 


40062 


3.92 


5.41 


10.61 


34.17 


990 


-6.152 


2762 


0.84 


yes 


-6.265 


3260 


0.77 


-6.011 


6 


38151 


3.92 


5.30 


10.23 


31.73 


994 


-6.386 


2697 


0.81 


yes 


-6.493 


3277 


0.77 


-6.234 


6.5 


36826 


3.92 


5.20 


9.79 


29.02 


983 


[-7.243] 


[6395] 


[1.87] 


no 


-6.918 


5244 


0.95 


-6.427 


7 


35531 


3.92 


5.10 


9.37 


26.52 


972 


[-7.340] 


[7325] 


[2.65] 


no 


- 


- 


- 


-6.624 


7.5 


34419 


3.92 


5.00 


8.94 


24.15 


959 


[-7.745] 


[12028] 


[1.81] 


no 


- 


- 


- 


-6.820 


8 


33383 


3.92 


4.90 


8.52 


21.95 


944 


[-7.781] 


[10650] 


[1.57] 


no 


- 


- 


- 


-7.019 


8.5 


32522 


3.92 


4.82 


8.11 


19.82 


923 


[-7.802] 


[9427] 


[1.40] 


no 


- 


- 


- 


-7.167 


9 


31524 


3.92 


4.72 


7.73 


18.03 


908 


[-7.818] 


[8283] 


[1.21] 


no 


— 


- 


- 


-7.374 


9.5 


30488 


3.92 


4.62 


7.39 


16.46 


892 


[-7.793] 


[6704] 


[1.10] 


no 








-7.590 


Giants 




























3 


42942 


3.77 


5.92 


16.57 


58.62 


915 


-5.445 


3275 


0.90 


yes 


-5.551 


3756 


0.87 


-5.182 


4 


41486 


3.73 


5.82 


15.83 


48.80 


866 


-5.540 


2945 


0.90 


yes 


-5.641 


3272 


0.84 


-5.303 


5 


39507 


3.69 


5.70 


15.26 


41.48 


837 


-5.630 


2460 


0.90 


yes 


-5.810 


3053 


0.83 


-5.491 


5.5 


38003 


3.67 


5.63 


15.13 


38.92 


833 


-5.867 


2852 


0.96 


yes 


-5.946 


3160 


0.83 


-5.629 


6 


36673 


3.65 


5.56 


14.97 


36.38 


825 


-6.100 


3165 


0.98 


yes 


-6.108 


3200 


0.84 


-5.769 


6.5 


35644 


3.63 


5.49 


14.74 


33.68 


810 


-6.278 


3534 


1.07 


yes 


-6.320 


3743 


0.91 


-5.902 


7 


34638 


3.61 


5.43 


14.51 


31.17 


798 


[-6.804] 


[7140] 


[3.46] 


no 


- 


- 


- 


-6.016 


7.5 


33487 


3.59 


5.36 


14.34 


29.06 


785 


-6.606 


4408 


1.20 


yes 


- 


- 


- 


-6.166 


8 


32573 


3.57 


5.30 


14.11 


26.89 


768 


-6.655 


3668 


1.05 


yes 


-6.692 


3857 


0.93 


-6.286 


8.5 


31689 


3.55 


5.24 


13.88 


24.84 


749 


-6.557 


2266 


0.80 


yes 


-6.770 


3490 


0.91 


-6.409 


9 


30737 


3.53 


5.17 


13.69 


23.07 


733 


-6.812 


2960 


0.90 


yes 




- 


- 


-6.564 


9.5 


30231 


3.51 


5.12 


13.37 


21.04 


709 


-6.848 


2594 


0.89 


yes 


-6.923 


3002 


0.85 


-6.646 


Supergiants 




























3 


42551 


3.73 


6.00 


18.47 


66.89 


912 


-5.347 


3346 


0.92 


yes 


-5.445 


3719 


0.86 


-5.083 


4 


40702 


3.65 


5.94 


18.91 


58.03 


837 


-5.387 


2877 


0.92 


yes 


-5.497 


3299 


0.86 


-5.144 


5 


38520 


3.57 


5.87 


19.48 


50.87 


779 


-5.561 


2974 


0.95 


yes 


-5.554 


3030 


0.86 


-5.247 


5.5 


37070 


3.52 


5.82 


19.92 


48.29 


764 


-5.611 


2938 


1.04 


yes 


-5.664 


3153 


0.87 


-5.352 


6 


35747 


3.48 


5.78 


20.33 


45.78 


747 


-5.751 


3000 


1.05 


yes 


-5.814 


3270 


0.90 


-5.438 


6.5 


34654 


3.44 


5.74 


20.68 


43.10 


732 


-5.945 


3531 


1.16 


yes 


-5.920 


3328 


0.93 


-5.520 


7 


33326 


3.40 


5.69 


21.14 


40.91 


715 


-5.995 


3230 


1.09 


yes 


-6.059 


3606 


0.96 


-5.642 


7.5 


31913 


3.36 


5.64 


21.69 


39.17 


702 


-6.036 


2702 


1.03 


yes 


-6.116 


3043 


0.90 


-5.781 


8 


31009 


3.32 


5.60 


22.03 


36.77 


678 


-6.058 


2366 


1.06 


yes 


-6.181 


2756 


0.88 


-5.873 


8.5 


30504 


3.28 


5.58 


22.20 


33.90 


644 


-6.143 


2498 


1.08 


yes 


-6.189 


2572 


0.90 


-5.895 


9 


29569 


3.23 


5.54 


22.60 


31.95 


629 


-6.385 


2988 


1.05 


yes 


-6.319 


2640 


0.91 


-5.998 


9.5 


28430 


3.19 


5.49 


23.11 


30.41 


613 


-6.487 


2921 


1.08 


yes 


-6.449 


2642 


0.93 


-6.148 



4.1. Comparison to Vink et al. mass-loss recipe 



The Monte Carlo method by Ide Koter et al. (Il997h as sum- 
marized in Sect. |2]has been used by Vink et al. (l2000l 1200 ih 
to compute a grid of mass-loss rates for O-type stars from 
which a fitting formula has been derived that provides M as 
a function of luminosity, effective temperature, mass and the 
ratio of the terminal velocity over the effective escape veloc- 
ity, i.e. Vca/vsic- This mass loss prescription is widely used in 
stellar evoluti on pre dictions (se e e.g iMevnet & Maeder i2003i: 
Palacios et al .l 120051: Em ongi & Chieffi[l2006l: lEldridge & VinkI 
'2006tiBrott et al Il2009l:1 Vink et al.ll201oh . 



For the spectral range that is investigated here the canonical 
value, based on empirical findings, for the ratio of t erminal ve- 
locity over escape velocity is 2.6. To compare to the I Vink et al.l 
( l2000h results, we calculated the mass-loss rate of our grid of 
stars using their prescription, that assumes [3-1. The results 
are given in the last column of Table [T] Figure |7] shows the to- 
tal energy that is extracted from the radiation field and that is 
transferred to th e stellar wi nd for all three methods: best-/?, hy- 
drodynamic and lVink et al.l prescription. All three methods yield 
similar, but not identical results in the regime where the best-/? 
and hydrodynamical method provide physical solutions. In terms 
of mass loss rates, we find that the predictions with the best-y6 
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Fig. 6. The Monte Carlo line force as a fraction of the effective 
gravity in the region around the sonic point for the best-/? so- 
lution of the 06 V (top panel) and 06.5 V (bottom panel) stars. 
The contribution of iron is shown separately (blue dotted line). 
Note that in the case of the 06.5 star the line acceleration does 
not balance the effective gravity at the sonic point. This is inter- 
preted as a failure to support a line driven wind. 



and hydro dynamical m ethod are on average about 0.2 to 0.3 dex 
lower than lVink et al.L again with the clear exception of the stars 
for which we fail to drive a stellar wind. As suggested by the 
similar wind energies the terminal velocities predicted by our 
best-yS and hydrodynam ical method turn out to be higher than 
adopted by IVink et aU We discus s these Voo in more detail in 
Sect. I4.3l as well as the reason whv lVink et al.l are able to predict 
M values for late O-type dwarfs and giants, where we fail. 

We emphasize that if the Wink et al] prescription is used as- 
suming the terminal velocities predicted by our best-f} or hydro- 
dynamical method, wherever these yield physical solutions, the 
mass loss rates agree to within ~0.1 dex. 



4.2. Comparison to (Modified)CAK-theory 



Since iLucv & SolomonI d 19701) it is generally accepted that the 
winds of massive stars are driven by the transfer of momentum 
(and energy) from the radiation field to the atmospheric gas, 
and that atomic tra nsitions play a pivotal role in this process. 
ICastor et al.l d 19751) describe the force associated to atomic tran- 
sitions by introducing a force multiplier 



e (f) 

M(^)-^-^^^ 

gcit) 



(19) 



where k and a are fitting parameters and t is an optical depth like 
parameter given by: 



/ du 

t ^ O-ePUthi — 



(20) 



Here p is the density, i>th the thermal v elocity of carbon ions 
at the effective temperature of the star dPauldrach et al.l [T986h 
and (Te the mass scattering coefficient of the free electrons. This 
parametrization of the line force is based on the expression of 
the force multiplier for a single spectral line. 



Mline(0=^^-[l-exp(77f)]: 

F t 



(21) 



where Avq is the Doppler shift of the frequency of the spectral 
line due to the thermal velocity of the particles in the wind, Fy is 
the flux at frequency v, F the total flux and rj is the ratio of line 
opacity to electron scattering opacity. Note that for optically thin 
lines Miine(0 becomes independent of f, whilst for optically thick 
lines Miine(f) <^ t. The cumulative effect of an ensemble of lines 
of various strengths is then expressed by Eq.[T9l The constant a 
in this expression is a measure of the ratio of line acceleration 
from optically thick lines only to the total line acceleration and 
k is related to the overall (line)strength of the ensemble of lines. 
See[Puls et al. (200(^ for a more in depth discussion on the CAK 
line force. It is a ssumed that a a nd k are constants throughout 
the wind (but see lKudritzkill2002h . The effects of changes in the 
ionization structure of the wind are modeled by multiplying ex - 
pression [T9lbv the term («e/W)'*, introduced bv lAbbottl (Il982h . 
«e is the electron number density, W the dilution factor and 6 is 
a constant. 

The parametrization of the line force as given by Eq.[T9lleads 
to an analytical expression for M and as a function of the fit- 
ting parameters k a nd a and the s t ellar p aramete rs. These e xpres- 
sions are given by ICastor et alJ (Il975i) . Pauldrach et al.l (Il986l) 
extend these expressions to account for the finite size of the stel- 
lar disk. 

To allow for a comparison with Abbott's results for the be- 
havior of the force multiplier as a function of f, we ignore the 
effect of n^/W. We calculated the force multiplier of the simu- 
lated line force, i.e. Eq. |6] at all our radius grid points and de- 
termined the corresponding value t. Typically, the optical depth 
like parameter ranges from t = lO"*" to large t. At large t the 
line force can be neglected com pared to the continuum radiation 
force. Following lAbbottl (119821) . we do not consider these large t 
points in this discussion but we focus on the range t < 10 

Figure [8] compares the behavior of the force multiplier for 
our best-j8 solution of the 03 V star Note that our Monte Carlo 
solution shows that M(t) can not be fully described by a strict 
power law, as assum ed in (mod ified) CAK, or equivalently, a is 
not inde pendent of f (IVinkl200(3|). Two causes ca n be pointed out 
(see e.g. |Puls"l987t ISchaerer & de Koteij[T997l) : /) the presence 
of a diffuse field due to multiple scatterings; ii) a complex behav- 
ior of radial stratification of the excitation and ionization, specif- 
ically near the sonic point. Here the wind accelerates rapidly, 
which causes a sudden steep drop in the electron density. As a 
result, elements which happen to have two dominant ionization 
stages (near r^) may temporarily re-ionize. 

The introduction of a ^-dependence of the CAK force multi- 
plier in Eq.[T9] by adding the term («e/ W)*, does not improve the 
fit to the Monte Carlo line force. This extended CAK description 
describes a plane through the three dimensional space spanned 
by the logarithms of t, n^/W and M(f), while the Monte Carlo 
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log T [kK] log T [kK] log T [kK] 

Fig. 7. The wind energy as a function of the effective temperature for dwarfs (left panel), giants (middle panel) and supergi ants (right 
panel). The best-/? method is given by the red squares, the hydrodynamic method by the green circles and the lVink et al.k ecipe by 
the blue triangles. Note that the kinetic energy in the wind is almost equal for all three methods. The mass-loss rate and terminal 
velocity for the three methods vary. 



line force follows a curved line through this space and is thus 
not confined to that plane. 

Our fitting function, Eq.|6l nicely captures the curved behav- 
ior of the line force in the supersonic part of the flow0 



- We also compared our fitting formula to the output of the starting 
model of the iteration cycle of the 03 V star. This allows us to assess 




10 12 14 16 18 20 
R[R.] 



Fig. 8. Top panel: The logarithm of the force multiplier M{t) for 
the hydrodynamic solution of our 03 V model (red plusses) as 
a function of the optical depth like parameter f; our fit function 
Eq. |6]to this data (green line) and the CAK fit function Eq. [19] 
to this data (blue line). Bottom panel: The force multiplier as a 
function of radius r for the same star Note that the CAK force 
multiplier is smaller than ours for large radii, resulting in a lower 
predicted Voo- 



Using the force multipliers k and a as derived from the 
Monte Carlo line force we can comp ute the CAK mass loss and 
terminal velocity dCastor et al .11 19751) . The Mcak values derived 
from these k and a are typically 0.0 to 0.3 dex higher than our 
best-/? results and our hydrodynamic solutions. A comparison of 
the terminal velocities is not meaningful since the slope in the su- 
personic part of the wind is not represented well by a. Therefore, 
the velocities derived with our k and a are on the order of the 
escape velocity. If w e compare to the mod ified CAK terminal 
velocities, following dPauldrach et al.|[T986h . we note that they 
are slightly lower than the velocities we derive (see also section 

4.3. Comparison with observations 

In this section, we compare our results to observations. We first 
compare predicted and empirical terminal velocities. Given that 
we find that for stars more luminous t han 10^-^ Lq our mass-loss 
rates agree well with the IVink et aTl prescription - whic h has 
been extensively scrutinized (see e.g.ljRepolust et al. 20041) - we 
focus the comparison of empirical and predicted mass-loss rates 
on lower luminosity stars, for which a 'weak-wind problem' has 
been identified. 



4.3.1. Terminal velocities 



been devoted to 
early- type stars. 



measuring the 
Summariz ing 



ter- 
the 



Several studies have 
minal vel ocities of 

work by |H owarth & Pri nial "(Il989h: iPriniaetalJ ([ 19901) : 
iLamers et airdl995.) : .Howarth et al.l (Il997h: iPuls et al.l (Il996l) . 
and iKudritzki et alj (1 19991) . iKudritzki & Pulsl (l2000h derive 
that the average value of empirically determined terminal 
velocities for stars hotter than 21000K is Uoo = 2.65 Uesc- The 
quoted accuracy of this mean value is roughly 20 percent. The 
Vco values are "measured" from the maximum blue-shifted 
absorption Umax in resonance Unes of ions such as Civ, Nv 
and Si iv, located in the ultraviolet part of the spectrum. These 
measurements are prone to systematic uncertainties, that have 
been extensive discussed in the literature (see for instance the 
above references). They may work in both directions. Effects 
that may cause the terminal velocity to be higher than Umax are 
measurements from lines that are not saturated in the outer wind 



whether the iteration procedure perhaps forces the line force into the 
shape of Eq. [6] However, also the line force resulting from the first 
iteration cycle is well represented by our description, indicating that it 
is quite generic. 
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4.55 4.6 
log T [K] 

Fig. 9. Predictions of the ratio of terminal velocity over effec- 
tive escape velocity at the surface for both the best-/? (in red) 
and hydrodynamic method (in green) as a functio n of effective 
tempe rature. The blue circles denote the data from lLamers et al.l 
(Il995h and the black tri angles show the empirical values from 
lHowarth&Prinial(ll989l) in which Vcx> is determined from the ul- 
traviolet P-Cygni profiles. 



(where for all practical purposes Vco is reached) or from ions 
that recombine in the outer wind. The former may be expected 
for stars with weak winds, the latter is more likely to occur 
in very dense winds. Effects that may cause Voc to be smaller 
than Umax may be the presence of turbulence in the outflow or 
the presence of strong atmospheric absorption at wavelengths 
slightly bluer than the wavelength corresponding to the terminal 
velocity, mistakenly contributed to absorption in the resonance 
line. Given the possible occurrence of these systematic effects, 
the uncertainty in the terminal velocity may be 10-15 percent for 
supergiants, and substantially larger than 20 percent for dwarfs. 
The error in the ratio yoo/fesc also includes uncertainties in (;e,sc- 
The largest contribution to this error comes from uncertainties 
in the stellar masses, that have been derived from a comparison 
to tracks of stellar evolution. It seems realistic to adopt a 30-40 
percent uncertainty in the empirical values of Uoa /v^sc rather than 
the 20 percent; quoted at the beginning of this section. 

Although iKudritzki & Puis! (l2000l) (and also iLamers et all 

[1991 conclud e that the ratio t^oo/ t^ esc is c onstant for O-type stars, 
the results of iHowarth & Priiiial (1 19891) show this ratio to de- 
crease with temperature, from about 3.5 at 31 500 K to about 2.4 
at 43 500 K. A luminosity class dependence of Hoo/fesc has to our 
knowledge not yet been reported. 

Figure |9] shows our predictions of Voc/vesc plotted against 
temperature. Different symbol types denote best-;6, hydrody- 
namic resu lts and emp i rical v alues from iHowarth & Prinial 
(119891) and iLamers et all (Il995h . As discussed in Sections |Z4| 
and 12. 5 1 our predictions of Uoa have individual random error bars 
of 10 to 20 percent 

In all cases, our predictions result in terminal velocities that 
are larger than observed. For the main sequence 03-06 stars the 
mean predicted ratio is 3.1 for the best-/? method and 3.6 for 
the hydrodynamical method. This is 17% and 36% higher than 
the observed mean value of 2.65. For giants the discrepancies 
are respectively 30% and 44%. For the supergiants the largest 
discrepancies are found. Using all O stars, the best-/? method 
over-predicts the ratio by 54%. The hydrodynamical method 
yields values that are on average 60% higher. The discrepancy 
between theory and observations thus seems to increase from 



dwarfs to giants to supergiants. Given the uncertainties, the over- 
prediction for the dwarfs may not be significant. 

The predictions show a tentative trend of a decreasing Uoc /fesc 
with temperature. As pointed out, recent empirical studies do not 
recover this behavior. Int erestingly, this type of tren d appears to 
be visible in the study bv IHowarth & Prinial (Il989h . Their trend 
is plotted in Fig. |9] featuring a slope that is comparable to the 
slope of our predictions. However, given the uncertainties in the 
current empirical estimates of Dm, we do not feel that this can be 
applied to (further) scrutinize the theory. 

Larger predicte d terminal velocities are also reported by 
ILamers et alJ (Il995h . In their sample, dominated by supergiants, 
a comparison to CAK models yields over-predictions by about 
33 percent, so slightly less compared to what we find. The rea- 
son for the over-predicted Vca values is unclear. Possible explana- 
tions (for part of the problem) include, /) overestimated correc- 
tions for the effect of turbulence (see above), a clumped and 
porous outer wind, hampering the acceleration of the flow in this 
part of the outflow from reaching as high a terminal velocity a s 
predicted here and in (modified) CAK (see lMuiires et al.ll201 lb . 
or ///) a systematic over-estimate of stellar masses. A system- 
atic discrepancy between masses of galactic stars derived from 
comparing their positions in the Hertzsprung-Russell diagram 
to evolutionary tracks and masses calculated fro m the spectro- 
scopically determined gravity was reported by e.g. lHerrero et al.l 
(1992) (but see Weidner & Vink 2010). Improvements in both 
the model atmospheres and fitting procedure seem to have re- 
duce d, but possibly not yet eliminated, the size of this discrep- 
ancy (iRepolust et al.ll2004l: iMokiem et al.ll2005l) . 

Unfortunately, progress in resolving the differences in pre- 
dicted and empirical yoo/uesc ratios quite strongly depends on our 
knowledge of stellar masses. Hopefully, detailed studies of very 
large pop ulations, such as th e VLT-FLAMES Survey of mas- 
sive sta rs (lEvans et al.ll2()05l) and the VLT-FLAMES Tarantula 
Survey ( lEvans et al. 1201 1|) may help resolve this issue. 



4.3.2. Mass loss rates: the weak-wind problem 

Relatively recent, analysis of appreciable samples of galactic 
stars using sophisticated model atmospheres has revealed a mis- 
match, possibly as high as a factor of 100, between empirically 
derived mass-loss rates and theoretical pr edictions for stars less 
luminous than abou t 10^^ Lp, (see e .g. Mart ins et al.l l2005bt 
iMokiem et all 120071: iMarcohno et"an 2009, and Fig. [TOll. As 
mass loss scales with some power of the luminosity, this problem 
occurs below a critical mass-flux an d is termed the 'weak- wind' 
problem (for a recent review, see iPuls et all |2008). Proposed 
explanations address deficiencies in determining the empirical 
mass-loss rates as well as in mass-loss predictions. Regarding 
empirical M determinations it should be realized that only UV 
resonance lines can be used as a diagnostic in the weak-wind 
regime, whilst in the (lets call it the) strong wind regime Ha and, 
in the Galactic case, radio-fluxes may also be used. The ions that 
produce the UV resonance profiles, such as Civ, N v and Siiv, 
often represent minor ionization species. The ionization continua 
of these species border the soft X-ray regime and therefore wind 
material may be susceptible to (non-thermal) processes produc- 
ing soft X-ray emissio n, such as shocks or magnetic mechanisms 
(iMartins et al.ll2005bl) . 

From a theoretical viewpoint, potential causes of the weak- 
wind problem include the decoupling of the major driving 
ions (the metals) from the bulk of the plasma at low densi- 
ties, whe n Coulomb coupling fails , and t he subsequen t ionic 
runaway (ISoringmann & PauldrachI [1991 iBabell I1995L 119961: 
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Krticka & Ktibi^ l2000t lOwocki & Pulsl l2002t iKrticka et all 
2003'): the shadowing of wind-driving lines by photospheric 
lines (Bab el 199 6j, and the neglect of curvature ter ms in the ve- 
locity field dPuls et al.lll998llOwocki & Pulslll999h . 

The results presented in this paper point to a cause for the 
weak winds related to the predictions of mass loss. This potential 
cause was quantitatively explored bv lLucvl (1201 Ol) . w ho pointed 
out tha t the global dynamical constraint imposed by IVink et all 
(I2OOOI) and recapped in Sec. 12.21 (notably Eq. HJi need not guar- 
antee that the derived mass-loss rates are consistent with station- 
ary trans-sonic flows. Here we h ave shown that although this 
assumption by IVink et all ( l2000l) is allowed for stars with lu- 
minosities above 10^'^ L©, it is not for lower luminosities. This 
luminosity limit for galactic O stars agrees with the empirical 
limit at ~ 10^'^ Lq to within 0. 1 dex in log L. The physical cause 
of the different M regimes (weak and strong winds) is a lack of 
line acceleration at the base of the wind. In main sequence stars, 
a contribution of Fe v lines is present in 06 stars and is miss- 
ing in the (lower luminosity and cooler) 06.5 stars, where Feiv 
is more dominant (see Sect. 13.3b . The importance of the Fe v/iv 
ionization balance has been pointed out by iLucvl (1201 Ol) and is 
confirmed by our results. 

In one f undamental aspect our results differ from that of 

ILucvI (l2010l) . For the sample o f low luminosity stars ( i. e. less 
than 1 0^-^ Lq) investigated by iMarcoUno et al.l (l2009i) . iLucvl 
(I2OIOI) predict s mass-loss r ates th at are about 1.4 dex lower than 
anticipated bv IVink et al.l (|2000() . The hydrodynamical method 
presented in this paper, which identifies the ability of the star to 
drive an outflow with a balance of the line force and the gravi- 
tational force at the sonic point, predicts that these stars do not 
have a wind at all. As it is clear form the presence of the shape of 
UV resonance lines that these stars do have stellar winds (with 
average mass loss rates that are 0.8 dex lower than iLucvl ' s pre- 
dictions), our result suggests that either some other mechanism 
is driving the wind or is supplementing the line acceleration at 
the base of the wind. Which force (or forces) counterbalances 
gravity remains to be identified, but perhaps magnetic pressure, 
effects of turbulence, and/or pulsations may play a role. 

We do point out that once material is accelerated, there 
is sufficient opacity avai lable to furth e r accele rate it to larger 
velocities. Interestingly, iMartins et al.l ( l2005bb report that for 
their sample of galactic weak-wind objects the average value of 
foo/t^esc is rather close to unity, and not 2.65. Although they con- 
cede that given the low wind densities their Voo values may be 
lower li mits, they po i nt to a mechanism of X-ray heating pro- 
posed by iDrew et all ( 1 19941) that may perhaps explain these re- 
sults. In the outer atmospheres of weak winds the cooling times 
can become quite long, such that heating of the material in for 
instance shocks may warm up the medium and strongly modify 
the ionization structure, in effect canceling the line force. 

Modified CAK theory does not predict the weak-wind dis- 
continuity. In this theory the adopted values for k and a are 
based on the input stellar spectrum, while the dilution and exci- 
tation/ionization changes throughout the wind are described by a 
fixed 6. Therefore, no self-consistent feedback between the wind 
properties and the li ne acceleration is accounted for. We note that 
in the predictions bv lPauldrach et al.l (1200 ll) a change of slope of 
the modified wind momentum luminosity relation can be seen at 
a luminosity of about 10^'^ Lq. It is tentative to suggest that if 
iPauldrach et al.l would have implemented the iterative procedure 
that we use, they might have identified a weak-wind regime on 
theoretical grounds. 




4.5 5.0 5.5 6.0 

log(L/Lsun) 

Fig. 10. The modified wind momentum D^om - ^R* IRp))Mvoo 
as a fu nction of stellar luminosity using the data of lMokiem et al.l 
( l2007h . Black symbols refer to mass-loss estimates based on the 
fitting of the Ho'-profile; grey symbols are mass-loss estimates 
that rely strongly on ultraviolet resonance lines. Note that the Ha 
estimates at L < IO^^Lq are upper limits. A steep jump-of about 
2 dex-can be seen at a luminosity of 10^ "L©. The red dots are 
our predictions for Dmom- The squares denote the supergiants, 
the circles the giants and the triangles the main-sequence stars. 
Below 10''^ Lq we do not find wind solutions for dwarf stars. At 
higher luminosity our predictions are close to the observed val- 
ues. We therefore interpret the origin of the weak-wind problem 
in dwarf stars to be connected to a lack of line driving for objects 
less bright than about 10^'^ . 

5. Conclusions 

We have presented new mass-loss rates and terminal wind veloc- 
ities for a grid of massive O-type stars, i mproving the treatmen t 
of p hysics in the Monte C arlo method bv I Abbott & Lucvl (Il985h 
and Ide Koter et al.l (1 19971) to predict wind properties of early- 
type stars. Two new type s of solutions have been discussed. First, 
building on the work of iMiiller & VinkI (l2008h . we present so- 
called best-yS solutions in which one still assumes a jS-type veloc- 
ity law (see Eq.|2]i in the wind, but in which the terminal veloc- 
ity and p are no longer adopted but constrained by requiring that 
they best fit the line force (distribution). Second, we abandon 
the ;6-type velocity structure and introduce numerical solutions 
of the wind stratification. Our main conclusions are; 

1. For stars more luminous than 10^'^ Lq the best-yS and hydro- 
dynamical method yield fi and Voo results in agreement with 
each other (within 5-20 percent), whilst the mass-loss rates 
agree within a factor of 2. 

2. Furthermore, both methods are in very good agree ment 
with the mass-loss prescription by IVink et aTl (l2000l) us- 
ing our terminal velocities in their recipe. This implies that 
the main assumption entering the method on which the 
IVink et al.l results are based - i.e. that the momentum equa- 
tion is not solved explicitly - is not compromising their pre- 
dicted M in this lumin osity range . Terminal velocity is an 
input parameter to the IVink et al] recipe. If we apply the 
canonical value Uoa = 2.6t;esc, the discrepancy between our 
mass-loss rates and their mass-loss rates is of the order of 
0.2, although occasionally 0.6 dex. 

3. At luminosities ^ 10^^ Lq our hydrodynamical method fails 
to produce an outflow because of a lack of line driving at the 
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base of the wind. This critical luminosity coincides with the 
onset of the 'weak- wind problem'. 

4. For O dwarfs the above luminosity criterion translates to 
a boundary between starting and failing to start a wind at 
06/06.5. The direct cause of the failure to start a wind in 
06.5 V stars is the lower luminosity and the lack of Fe v 
lines at the base of the wind compared to 06 V stars. 

5. The fact that our hydrodynamical method fails to drive a 
wind at L ^ 10^'^ Lq may imply that some other mechanism 
is driving the weak winds or is supplementing the line accel- 
eration at the base of the wind to help drive gas and initiate 
the wind. 

6. For stars more luminous than 10^^ we predict, using the 
best-yS and hydrodynamical method, terminal velocities that 
are typically 35 and 45 percent higher than observed. Such 
over-predictions are similar to what is seen in MCAK-theory 
jLamers et al.lll995h . 

7. We predict beta values in the range 0.85 to 1 .05, with a trend 
that supergiants have slightly higher /3 values than dwarfs. 
Th is range of P values agrees very well with empirical results 
bv lMassaet all (120031) . 
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